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INTRAOPERATIVE MONITORING OF TEMPERATURE-INDUCED 
TISSUE CHANGES WITH A HIGH-RESOLUTION 
DIGITAL X-RAY SYSTEM DURING THERMOTHERAPY 



The U.S. Government has a paid-up license in this invention and the right in limited 
10 circumstances to require the patent owner to license others on reasonable terms as provided for 
by the terms of Grant or Contract No.: USAMRDC/AIBS#2459, log #B4340289, 9/15/94 to 
9/14/96, awarded by the Department of Defense. 

RELATED APPLICATION 
This application is a continuation of U.S. Patent Application No. 09/556,958 filed on 
1 5 April 2 1 , 2000, the entirety of which is hereby incorporated by reference. 

BACKGROUND OF THE INVENTION 
L Field of the Invention 

This invention pertains to intraoperative monitoring methods and devices, and more 
particularly to intraoperative monitoring devices and methods for monitoring temperature- 
20 induced changes to tissue with high resolution digital X-ray imaging and inducing the changes to 
the tissue. 

2 Discussion of Related Art 

Various types of thermotherapy have been considered, and/or applied to the treatment of 
cancers. Laser interstitial thermotherapy (LITT) has received attention as a surgical procedure 
25 for the treatment of cancer tumors, and particularly with respect to liver, head, neck and breast 
cancer. 

Laser interstitial thermotherapy is a surgical procedure for the treatment of cancer tumors 
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5 in which near infrared laser energy is delivered to the tumor site inside the body through a 

flexible fiber-optic probe. Some fiber optic probes in use terminate with a light diffusing tip. The 
infra-red laser radiation is absorbed by the tumor cells, which results in a temperature increase 
and subsequent cell death of the tumor cells. The temperature distribution around the light- 
diffusing tip, and thus the extent of cell-death, is a function of the laser parameters, treatment 

10 time, tumor size, and shape, fiber optic tip geometry, optical properties of the tumor, and blood 
perfusion rates in both compressed and uncompressed treatment sites. 

The treatment parameters (e.g., wavelength, power, duration, tip geometry, and tip 
location and orientation) must be selected so as to minimize collateral damage to healthy tissue 
surrounding the tumor, yet still must ensure reliable total tumor destruction. Because of tissue 

1 5 inhomogeneities, and inter-patient variability of the physical and biological properties of tumors, 
intraoperative monitoring of the treatment effect is highly desirable. Currently, intraoperative . 
monitoring of LITT is conducted with magnetic resonance imaging (MRT) or three-dimensional 
ultra-sound, or by measuring the temperature at discrete locations in-situ with thermo-couples or 
thermo-sensing fluorescent probes. 

20 Intraoperative monitoring with MRI has numerous disadvantages which include being 

expensive and not being able to be used with metal protected light guides which are currently 
used for the LITT probe. Three-dimensional ultra-sound imaging techniques are currently at an 
experimental stage and have not been sufficiently developed. In situ thermo-couples or 
fluorescence-based temperature probes only provide temperatures at a relatively small number of 

25 points throughout the tissue being monitored. Currently, only fluorescence-based temperature 
probes are approved by the FDA for clinical use. 
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SUMMARY OF THE INVENTION 



It is thus an object of this invention to provide an intraoperative monitoring device and 



method used in conjunction with diagnostic procedures. 



It is another object of this invention to provide an intraoperative monitoring device and 



method used in conjunction with stereotactic X-ray mammography. 
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It is another object of this invention to provide an intraoperative monitoring method and » 



device that provides a high resolution X-ray image of temperature-induced changes to tissue 
during thermotherapy. 

It is another object of this invention to provide an intraoperative monitoring device and 
method that provides a high resolution temperature map of tissue during thermotherapy. ' •■ : \ 



real-time monitoring of temperature induced changes to tissue during thermotherapy using X-ray 

imaging for feedback information used during the thermotherapy. 

The above, and related objects of this invention are realized by providing a device for 

monitoring thermally-induced changes to localized regions of tissue which has an X-ray 
20 illumination source, an X-ray detector, a data storage unit in communication with the X-ray -i 

detector, an image comparison unit in comparison with at least the data storage unit and an ; 

image display unit in communication with the image comparison unit. The X-ray illumination 

source and the X-ray detector are arranged to reserve a space therebetween for. accommodating 

tissue to be monitored. Preferably, the tissue to be monitored is a portion of a patient's body 
25 which is being monitored during the surgical procedure. The X-ray detector produces a plurality 

of X-ray image data signals, each of which corresponds to an X-ray image of the portion of the 
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It is another object of this invention to provide a method of thermotherapy which includes 
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5 patient's body being monitored. Preferably, the X-ray detector produces a digital X-ray image 
data signal. The X-ray image data signal corresponds to a two-dimensional image of the portion 
of the patient's body being monitored in a preferred embodiment, and corresponds to a three- 
dimensional X-ray image of the portion of the patient's body in another preferred embodiment. 
The image comparison unit compares X-ray image values between corresponding spatial 
10 points, of first and second X-ray image signals to provide a resultant X-ray image signal based on 
the comparison. Preferably, the image comparison unit subtracts pixel values between 
corresponding spatial points of first and second X-ray image signals, providing a measure of the 
>.: : change in intensity of the received X-ray signal at each point within the digital 'X-ray image 
*: a. signals. The resultant image signal is then one particular example of a difference image signal 
15 that is generated by the image comparison unit and then displayed on an image display unit to 
; . provide real-time information concerning the temperature distribution and changes in 
temperature throughout the portion of the patient's body being monitored. The images also 
provide information corresponding to the volume of denatured tissue. In the preferred 
embodiment, both the data storage unit and image comparison unit are implemented within a 
20 i personal computer or workstation. * < • ? 

Another preferred embodiment of the present invention is directed to a device for causing 
thermally-induced changes to localized regions of tissue. The device according to this preferred 
> embodiment has an X-ray illumination source, an X-ray detector, a data storage, unit in 
communication with the X-ray detector, a thermotherapy heating assembly, an image 
25 comparison unit in communication with at least the data storage unit and an image display unit in 
communication with the image comparison unit. The combination of X-ray illumination source, 
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5 X-ray detector, data storage unit, image comparison unit, and image display unit are constructed 
and arranged in a manner similar to the monitoring device summarized above. The 
thermotherapy heating assembly may be selected from currently known devices and may include 
a laser irradiation devices, microwave irradiation devices, radio frequency irradiation device, or 
an ultra-sound energy source. In the preferred embodiment, the thermotherapy heating assembly 

L0 is a laser interstitial thermotherapy assembly. There are laser interstitial thermotherapy 
assemblies known in the art that are suitable for use with the device for causing thermally- 
induced changes to localized regions of tissue, in accordance with this invention. For example, 
the laser interstitial thermotherapy devices described in Robinson, David S. et al, "Interstitial 
Laser Hyperthermia Model Development for Minimally Invasive Therapy of Breast Carcinoma", 

15 J. Am Coil Surg, 1998, reprint pages 284-292; Milne, Peter J. et al, "Development of 

Stereotactically Guided Laser Interstitial Thermotherapy of Breast Cancer: In Situ Measurement 
and Analysis of the Temperature Field in Ex Vivo and In Vivo Adipose Tissue," 1 5 Lasers in 
Surgery and Medicine, 2000, reprint 26:67-75; and Maims, Fabrice et al, "In Situ Temperature 
Measurements With Thermocouple Probes During Laser Interstitial Thermotherapy (LTT): 

20 Quantification and Correction of a measurement Artifact," Lasers in Surgery and Medicine, Vol. 
23, No. 2, 1998, reprint pages 94-103 are suitable: the entire contentof each is incorporated 
herein by references. 

Another preferred embodiment of this invention is directed to a method of thermally 
inducing and monitoring changes to localized regions of tissue, including illuminating a volume 
25 of tissue with a first beam of X-rays, detecting portions of the first beam of X-rays that pass 
through localized regions of tissue within the volume of tissue, generating a first X-ray image 
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5 signal from the portions of the first X-ray beam detected, and applying heat to at least a localized 
region of tissue within the volume of tissue. A preferred embodiment of the method is directed 
to thermally inducing and monitoring changes to tumors in breast cancer patients. After . 
applying the heat, the volume of tissue is illuminated with a second beam of X-rays. X-rays from 
a second beam of X-rays that pass through the volume of tissue that includes the localized 
10 regions, e.g., through the tumors, are detected and a second X-ray image signal is generated 
therefrom. In the preferred embodiment, the first X-ray image signal is stored in a data storage 
unit and then retrieved for generating a resultant image signal which preferably is a d 
image signal based upon a comparison of the first and second X-ray image signals. In a 
■ J - preferred embodiment, the first and second X-ray image signals and the. difference image signal 
15 * are digital signals and the difference image signal is formed by subtracting each pixel value of 
the first X-ray image signal from a spatially corresponding pixel of the second X-ray image 
signal/ Each X-ray image signal is correlated with the detected X-ray intensity corresponding to 
: ■ the particular spatial point, and the difference image signal corresponds to an intensity change in 
-V the detected X-rays for each corresponding difference image point. The difference image signal 
20 is rendered as a difference image and displayed on an image display device, preferably; to i 
provide real-time feedback to the surgeon applying thermotherapy. The surgeon can then 
determine whether to alter the thermotherapy parameters, maintain the thermotherapy, or 
terminate the thermotherapy with the aid of the temperature change information displayed on the 
image display device. 

25 If the surgeon decides to continue the thermotherapy, the volume of tissue is illuminated 

with a third beam of X-rays. (The X-ray beams may be along the same or modified paths 
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5 relative to the volume of tissue.) The X-rays from the third beam that pass through the volume 
of tissue are detected and a third X-ray image signal is generated therefrom: The second and 
third beam of X-rays may be separated in time by a period in which there is no illumination of 
X-rays, or it may be a continuous illumination classified as contiguous time periods. Continuous 
illumination is currently less preferable than intermittent illumination due to safety concerns 
10 regarding the total X-ray dose applied to the patient. 

The first X-ray image signal is again retrieved from the data storage unit and used as a 
static reference image signal to produce a second difference image signal by subtracting each 
corresponding pixel of the first and second image signals. The second difference image signal is 
then rendered and displayed as an updated difference image which provides updated information 
1 5 to the surgeon, preferably in real time. The surgeon can then use the temperature change 

information displayed to reassess the status of the thermotherapy to determine whether to alter, 
continue or terminate the thermotherapy. 

This process is repeated until the surgeon determines that the thermotherapy should be 
terminated. In this embodiment, the difference image signal is always generated by retrieving the 
20 same static reference image signal from the data storage unit and subtracting it from the updated 
measurement signal. 

The above detailed description of a succession of measurements is by way of example. 
The reader should recognize from the teachings herein that the scope and spirit of the invention 
includes the general concepts and not the particular order of observing and responding to the 
25 observations. 

In an alternative embodiment, all method steps are the same as those 20 noted above, 
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5 except that a static reference image signal is not used to generate the difference image signal. In 
this preferred embodiment, the first X-ray image signal is replaced by the second X-ray image 
signal after the first difference image signal is generated. Similarly, the second X-ray image 
signal stored in the data storage unit is replaced with the third X-ray image signal after the 
second difference image signal is generated. This process is repeated until the surgeon 

10 determines that the thermotherapy should be terminated. This case provides a dynamic reference 
image signal for forming the difference image signal in which the dynamic reference image 
signal is updated after each succeeding illumination. 

In alternative embodiments, one may combine both static and dynamic processes for 
generating the difference image signals. For example, the reference image signal may be 

1 5 maintained in memory without being replaced for one, two or more successive illuminations, 

followed by being updated either frequently, such as with a pure dynamic reference image signal, 
or intermittently, again being a mix of dynamic and static processes. 

Another embodiment of the invention is directed to a method of destroying cancerous 
tissue by forming a first X-ray image of a portion of a patient's body, applying heat to a localized 

20 region of the portion of the patient's body, forming a second X-ray image of the portion of the 
patient's body subsequent to applying heat to the localized region of the portion of the patient's 
body, and generating a difference image based on a comparison of the first X-ray image data to 
the second X-ray image data. The surgeon then modifies the application of heat based on 
information obtained from the difference image. In a preferred embodiment, the first and second 

25 X-ray images and the difference X-ray image are high resolution, three-dimensional digital X- 
ray images. Preferably, the comparison is a subtraction of the first X-ray image from the second 
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5 X-ray image. The illumination with successive X-ray beams, detecting, generating X-ray image 
signals, and generating successive difference image signals is repeated as the surgeon requires 
until he terminates the thermotherapy. 

BRIEF DESCRIPTION OF THE DRAWINGS 

10 These and other objects and advantages of the invention will become more apparent and 

more readily appreciated, from the following detailed description of the presently preferred 
exemplary embodiments of the invention taken in conjunction with the accompanying drawings, 
of which: 

FIGURE 1 is a schematic illustration for explaining general concepts of the invention; 
1 5 FIGURE 2 is a schematic illustration of a device for causing and monitoring thermally 

induced changes to localized regions of tissue in accordance with a preferred embodiment of the 
invention; 

FIGURE 3 A is a partially cut-away view of a probe within a heating 5 assembly in a 
preferred embodiment of this invention; 
20 FIGURE 3B is a side view of an end of the probe illustrated, in FIGURE 3A; 

FIGURE 4 is a schematic illustration of applying heat to a tumor in accordance with the 
preferred embodiment of this invention; 

FIGURE 5 is a flowchart illustrating a preferred embodiment of a method according to 
this invention; 

25 FIGURE 6 is a flowchart illustrating another preferred embodiment of a method 

according to the invention; 
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FIGURE 7 is a graph showing average X-ray intensity in the region of 15 interest 

(diffusing fiber experiment); 

FIGURE 8 is a graph showing average X-ray intensity in the region of interest (bare fiber 

experiment); 

FIGURE 9 is a schematic illustration of the region of interest with respect to a 
thermotherapy probe for the experiments of FIGURES 7 and 8; 

FIGURE 10 shows graphs of changes of average X-ray intensity as a function of time for 
experimental regions of interest as illustrated in FIGURE 9; . ■ * 

FIGURE 1 1 illustrates desired hypothermic effects in accordance with the invention; 

FIGURE 12 is a plot of intensity data as a function of time over entire images taken every 
thirty seconds during an experiment in accordance with the invention; 

FIGURE 13 is a plot of intensity data as a function of time over whole 5 images taken 
every six seconds in accordance with the invention; 

FIGURE 14 is a plot of maximum pixel intensities in an experiment in accordance with 
the invention for thirty second firing sequence of the laser; 

FIGURE 15 is a graph of maximum pixel intensities as a function of time in a second 
experiment for six second firing sequences of the laser; 

FIGURE 16A illustrates transfer function as a function of intensity on a linear scale; 

FIGURE 16B illustrates a transfer function vs. intensity change on a logarithmic scale; 

FIGURE 16C illustrates a transfer function vs. intensity change for a 15 desired function; 

FIGURE 17A is an initial subtracted image on a linear scale for data taking in accordance 
with a preferred embodiment of this invention; 
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FIGURE 17B is an initial subtracted image on a logarithmic scale for data taken in 
accordance with a preferred embodiment of this invention; 

FIGURE 18A is a subtracted image after 3 minutes of laser firing on the linear scale for 
data taken in accordance with a preferred embodiment of this invention (difference image); 

FIGURE 18B is a subtracted image on a logarithmic scale corresponding to the case of 
FIGURE 18A (difference image); 

FIGURE 1 9 A is a subtracted image after 7 minutes of laser firing on a linear scale for 
data taken in accordance with a preferred embodiment of this invention (difference image); • . ■ . 

FIGURE 1 9B is a subtracted image after 7 minutes of laser firing on a 5 logarithmic : 
scale for the case corresponding to FIGURE 19A (difference 
image); 

FIGURE 20 A is a subtracted image after 4 minutes of cooling down on a linear scale for 
data taken in accordance with a preferred embodiment of this invention (difference image); 

FIGURE 20B is a subtracted image on a logarithmic scale for the case corresponding to 
FIGURE 20A (difference image); 

FIGURE 21 A is a subtracted image after 7 minutes of cooling down on a linear scale for 
data taken in accordance with a preferred embodiment of this invention (difference image); and 

FIGURE 21B is a subtracted image, on a logarithmic scale corresponding to the case of^ 
FIGURE 21 A (difference image). 
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5 DETAILED DESCRIPTION OF THE PREFERRED 

EMBODIMENTS 

We first explain the general concepts of the invention with reference to 20 FIGURE 1 . 
FIGURE 1 is a schematic illustration of a volume of tissue 100 which contains a localized region 

10 of tissue 102. In the preferred embodiments, the localized region of tissue 102 is a tumor formed 
of cancerous tissue. The volume of tissue 100 is illuminated with a collimated beam of X-rays 
104 illustrated schematically in FIGURE 1 as a heavy arrow and also labeled Io representing the 
intensity of the collimated X-ray beam. The collimated X-ray beam 104 has a plurality of beam 
portions, one of which is labeled 106 in FIGURE 1, that pass through a portion of the volume of 

15 tissue 100 for a distance 108 before reaching the depth of the tumor 102. The tumor 102 is 

represented as having a simple structure in FIGURE 1 for the purpose of illustrating the general 
concept of the invention. For facilitating the explanation, the tumor 102 has a uniform thickness 
1 1 0 along the direction of the path of the X-ray beam 1 04. In FIGURE 1 , the X-ray beam is 
illustrated as a collimated beam to facilitate the discussion. The X-ray beams emitted by sources 

20 currently used are diverging beams. Portions of the X-ray beam such as 106 pass through the 
volume of tissue 100 without passing through the lOtumor 102 and emerge from the volume of 
tissue 100 at an opposing side 112 of the volume of tissue 100. Again, for facilitating the 
explanation of the general concepts of the invention, the volume of tissue 100 is also assumed to 
have a simple structure in which it has a uniform width which is the sum of the distances 108 to 

25 reach the tumor 102, plus the distance 1 10 in order to pass by or through the tumor 102, plus the 
distance 1 14 to pass from an opposing end 1 16 of the tumor 102 to exit the volume of tissue 100 
at the opposing end 1 12 of the volume of tissue 100. 

It is well established that electromagnetic radiation has a dual nature in which it can be 
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5 thought of as a collection of a plurality of photons, or may be thought of as an electromagnetic 
wave or waves. The particle view is convenient to imagine that the X-ray beam is a collection of 
a large number of particles, i.e., photons, in which some of the plurality of photons traverse a 
path which misses the tumor 102, while another plurality of photons traverse a path through the 
volume of tissue 100 and then pass through the tumor 102. 25 Photons traveling along a path 106 

10 miss the tumor 102 and exit the volume of tissue 100 at the opposing end 1 12 after traveling 
through a distance of tissue which is the sum of the distances 108, plus 1 10, plus 1 14. In this 
illustrative example, the collection of all photons which pass above the tumor 102 in FIGURE 1 
and exit the opposing end 1 12 of the volume of tissue 100 together form an intensity sub-beam 
1 18, as illustrated in FIGURE 1, without passing through a tumor 102. Similarly, a portion of the 

1 5 X-ray beam 104 that passes through the volume of tissue 100 to emerge at the opposing end 112 
below the tumor 102, as illustrated in FIGURE 1, without passing through the tumor 102 forms a 
sub-beam intensity 120. A portion of the X-ray beam 104 which passes through a length of tissue 
108 before reaching the tumor 102, and passes through the length 1 10 of tumor 102, followed by 
passing through another length of tissue 1 14, emerges as an X-ray beam component 122. 

20 Conventional X-ray observation techniques rely on the fact that the tumor tissue 102 

attenuates and/or absorbs X-rays more than the surrounding tissue in the volume of tissue 100. 
Consequently, if the X-ray beam 104 entering the volume 100 is uniform when it enters the 
volume, the beam portions 118 and 120 will be more intense than the beam portion 122, due to 
additional attenuation and/or absorption by the tumor 102. The inventors recognized that the 

25 temperature-dependent attenuation and/or absorption of X-rays by tissue can be utilized to 

provide improved diagnostic methods and devices as well as to provide methods and devices for 
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5 a combination of therapeutic and diagnostic techniques. 

The above-described schematic illustration can also be modeled mathematically. The 
incident intensity is represented as Io and the detected intensity at a position x, y, in the plane 
perpendicular to the figure at opposing end 1 12 of the volume of tissue 100, is represented as 

i d {x>y) = he' lMw)dz (1) 

10 where ju (x,y,z) is the attenuation coefficient at each three-dimensional point within the volume 
of tissue. 

Upon selecting a coordinate system in which the plane orthogonal to the plane of 
FIGURE 1 at the opposing end 1 12 of the volume of tissue 100 is an X-Y plane and the Z- 
direction is the direction of travel of the X-ray beam, the lengths 108, 110 and 1 14 are 
15 represented as z\ 9 zz and Z3, respectively. In the model illustrated in FIGURE 1, a photon passing 
along the line 106 and exiting the volume of tissue 100 at the opposing end 1 12 is assumed to 
traverse a substantially homogeneous region of tissue in which the attenuation coefficient is 
constant in the Z-direction, and represented by ~ The intensity of X-rays passing through a 
portion of the volume of tissue 100 within the region 1 18 is then represented as 

20 ' I 2 =I 0 e-^ +Z ^ (2) 

where the integration performed along the Z-direction leads to the factor Zi+z 2 +z 3 in the 
exponential which is simply the total distance traveled through the volume of tissue 100 in the Z- 
direction. Similarly, the intensity of X-rays emerging at the region 122 that pass through two 
regions of normal tissue, as well as passing through the tumor, is represented as 
25 /i=V -^(W-^ 2 (3) 
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Where jj. 2 is the attenuation coefficient within the tumor which is substantially constant along 
the Z-direction, the attenuation coefficient ju 2 is generally different from the attenuation 
coefficient ju \. 

Generally, the tumor 102 attenuates X-rays more than that of surrounding healthy tissue. 
10 This is the case when the tumor is breast cancer. In that case, to a first approximation, one can 
assume that the healthy tissue causes essentially no attenuation to the beam of X-rays, while X- 
rays passing through the tumor are attenuated more strongly. In this approximation, equations 2 
and 3 become 

Wo (4) 

15 I } =I 0 e~ M ° Z2 e~ bTz2 (5) 

respectively. 

The inventors recognized that certain types of tissue have a temperature variation in its 
attenuation of X-rays. This can be represented in the above-noted illustrative model by an 
attenuation coefficient which varies with temperature. For example, a linear variation with 
20 temperature of the temperature coefficient ju = // 0 + bl leads to a temperature-dependent 
intensity 

7 1 (T) = V" // ° Z2 e" 6Tz2 (6) 

25 at the region 122. As one can see from equation 6, the intensity I- varies with temperature. If 
there was no variation in X-ray attenuation properties of the tumor with respect to changes in 
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5 temperature, then b would be equal to zero, and the second factor in equation 6 would be equal to 

unity. In that case, comparing the intensity of X-rays sent through the tumor at different 

temperatures would lead to substantially the same intensity. However, since there is a 

temperature variation in the attenuation properties of X-rays by the tumor, the intensity 

according to equation 6 at a temperature Ti is different from that of another temperature T 2 . This 

10 temperature-dependent variation in X-ray attenuation by the tumor can be exploited by observing 

the change in intensity I- as the temperature changes. For example, the ratio of the intensity at 

temperature T 2 to the intensity at temperature Ti, according to equation 6 can be represented as 

/ L (T 1 ) = ^z 2 (T 2 -T l ) 
/.fli) 

If there were no temperature variation in the X-ray attenuation coefficient of the tumor 
15 102, equation 7 would be equal to unity. The deviations from unity signify the temperature 

variations in the intensity of X-rays passing through the tumor 102. The ratio of intensities in this 
case is taken as a measure to recognize the change in intensity of X-rays passing through the 
tumor 102 with changes in temperature. However, the invention includes generally comparing 
; the changes in intensity of X-rays passing through tumors that are correlated with changes in 
20 k temperature and changes induced by temperature. For example, another useful method of 
comparing the changes of intensity of X-rays passing through the tumor with changes in 
temperature is to subtract the intensity values. There are innumerable measures that one may use 
to implement the general concept of comparing the intensity of X-rays passing through the tumor 
at one temperature to the intensity of X-rays passing through the tumor at another temperature. 
25 This invention is directed generally to the concept of utilizing such a change in X-ray intensity 
with a change in temperature in methods and devices. 
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5 The following detailed description of the preferred embodiments of the methods and 

apparatuses of this invention will refer particularly to the application of laser interstitial 
thermotherapy (LITT) which is a surgical procedure for the treatment of cancer tumors where 
near infra-red laser energy is delivered to the tumor site inside the body through a flexible fiber- 
optic probe that has a light diffusing tip. The inventors have found that the temperature increase 

10 generated during LITT will cause a change in the X-ray density of the heated tissue which is 
. detectable in accordance with this invention. In addition, thermal denaturation causes a variation 
in tissue X-ray density which is detectable according to this invention. The LITT procedure will 
be described in detail in the preferred embodiment; however, the general concepts of the 
invention include other mechanisms now known, and later developed which lead to heating 

15 tumor tissue. For example, lasers which irradiate tissue with electromagnetic radiation in regions 
other than the near infra-red region, microwave sources of radiation, radio frequency sources, 
and ultra-sound energy sources are all mechanisms that can be used to heat tumor tissue. 
Furthermore, the heat sources may be internally applied heat sources such as the LITT method, 
or they may be external sources. 

20 In addition, the preferred embodiments describe applications to breast cancer in particular 

detail. However, it is anticipated that other forms of tissue monitoring and treatment is included 
within the general concept of the invention. For example, the inventors anticipate that choroidal 
tumors of the eye, prostate cancer, and liver cancer are also particularly suitable for tumor 
monitoring and treatment according to this invention. The inventors anticipate that this invention 

25 is generally applicable to monitoring temperature changes in localized regions of tissue 
whenever the localized regions of tissue have a temperature variation in X-ray attenuation 
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5 compared to a surrounding volume of tissue. 

An example of a thermal therapy device in accordance with a preferred embodiment of 
this invention is illustrated in FIGURE 2, and labeled generally by reference numeral 124. The 
thermal therapy device 124 has a monitoring unit 126 and a heating assembly 128. The 
monitoring unit 126 has an X-ray illumination source 130 and an X-ray detector 132 disposed 

10 proximate to the X-ray illumination source 130. A space is reserved between an X-ray 

illumination source 130 and X-ray detector 132 for accommodating tissue to be monitored, 
and/or treated, such as a portion of a person's body. In the example illustrated in FIGURE 2, a 
woman's breast 134 is arranged between the X-ray illumination source 130 and the X-ray 
detector 132 for the treatment and treatment monitoring of breast cancer. A compression plate 

15 136 compresses the patient's breast 134 in contact with the X-ray detector 132. The monitoring 
unit 126 has a data storage unit 138 that is in communication with the X-ray detector 132. The 
monitoring unit 126 also has an image comparison unit 140 that is in communication with at 
least the data storage unit 138 (see, the dashed line in FIGURE 2). In an embodiment of the 
invention, the image comparison unit 140 is in communication with both the X-ray detector 132 

20 and the data storage unit 138. The image comparison unit 140 is in communication with an 

image display unit 142. The image display unit 142 may be a separate specially designed display 
unit, or may be the monitor of a personal computer or workstation 144. Each of the image 
comparison unit 140 and data storage unit 138 maybe separate dedicated components, or may be 
implemented on the personal computer or workstation. Furthermore, the image comparison unit 

25 140 maybe a stand-alone special function component, or may be implemented by programming 
the personal computer 144. Special function stand-alone units typically allow one to optimize 
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5 performance, but at a higher cost. 

In the preferred embodiment, the X-ray detector 132 is a digital detector that produces a 
digital X-ray image signal. The invention is not limited to digital X-ray detectors, and includes 
analog detectors producing analog X-ray image signals. In the preferred embodiment, the X-ray 
detector 132 has a fluorescent screen 146 which fluoresces in the visible region of the 

10 electromagnetic spectrum upon being struck by X-rays from X-ray illumination source 130. 

The light emitted by the fluorescent screen is detected by an array of CCD elements in 
the preferred embodiment. Furthermore, a preferred embodiment will typically have optical 
components arranged between the fluorescent screen and the CCD array (not shown in the 
drawings) in order to condense and focus the fluorescing light onto the CCD detector. There are 

15 high-resolution digital X-ray imaging devices currently on the market for treating and monitoring 
breast cancer which are suitable for use according to this invention. Preferably, the monitoring 
unit 126 employs a three-dimensional X-ray imaging system, such as stereotactic X-ray imaging. 
The MAMMO VISION- system 85200G-2 produced by Fischer Imaging Company was found to 
be suitable for use in the monitoring unit 126. The U.S. Patents Nos. 5,078,142; 5,365,562; 

20 5,526,394; and 5,9 1 7,88 1 , each of which is incorporated by reference herein in their entirety, 
describe details of various aspects of X-ray imaging systems suitable for application in the 
monitoring unit 126. 

In the preferred embodiment, the heating assembly 128 is a laser interstitial thermal 
therapy assembly. The LITT assembly 128 has a MAMMOGRAPHIC thermal ablative therapy 
- 25 probe 148 attached to sled 150 which is controlled by a data acquisition system and control unit 
152. An optical fiber 154 attached to the probe 148 to transport the laser energy from the laser 
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5 156. In a preferred embodiment, the laser is a continuous wave Nd:YAG laser. The invention 
anticipates that other lasers such as semiconductor lasers and gas lasers will be used in particular 
applications. Furthermore, it is anticipated that pulsed lasers may be suitable in some 
applications. The LITT assembly 128 also has a plurality of thermocouples 158. The insertion 
portion of the probe 148 is shown in more detail in FIGURES 3 A and 3B. The outer layers of the 

10 probe 148 are cut-away in cross-section in FIGURE 3 A to show a cross-sectional view of the 
quartz diffusing cap 160 and a cross-sectional view of the stainless steel tube 162. The probe 
148 has two thermocouples 164 and 166, respectively. The optical fibre 168 is exposed in the 
cut-away view of the quartz diffusing cap 160 and stainless steel tube 162. FIGURE 3B shows a 
side view corresponding to FIGURE 3 A. FIGURE 4 is a schematic illustration of the probe 148 

15 v inserted into a tumor 170 within a surrounding volume of tissue 172. 

FIGURE 5 illustrates a preferred embodiment of a method of thermally inducing and 
monitoring changes to localized regions of tissue, in accordance with the invention. A region of 
tissue 134 arranged between the X-ray illumination source 130 and the X-ray detector 142 is 
illuminated with a first beam of X-rays (see, FIGURE 2 in conjunction with FIGURE 5). The X- 

20 ray detector 1 32 detects a plurality of portions of the first beam of X-rays after they have passed 
through the volume of tissue 134. The X-ray detector 132 generates a first X-ray image signal 
from a plurality of portions of X-rays detected from the first beam of X-rays. In a preferred 
embodiment, the first X-ray image signal is a digital signal. However, the general concept of the 
invention includes generating an analog X-ray image signal, rather than a digital X-ray image 

25 signal. 

In a first preferred embodiment of a method of thermally inducing and monitoring 
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5 changes to localized regions of tissue in accordance with the invention, the first X-ray image 
signal generated is a static reference X-ray image signal which is stored in a data storage unit 
138 for later retrieval. 

Heat is then applied to the tissue using an appropriate heating assembly for thermal 
therapy. For example, heat may be applied with the LITT assembly 128 illustrated in FIGURE 2. 
10 However, as noted above, the heat application assembly may be selected according to the 
particular, use from known heating assemblies such as various laser assemblies, microwave 
heating, radio frequency heating or ultrasound heating assemblies. 

The tissue 134 is then illuminated with a second beam of X-rays at a second time. The X- 
ray detector 132 detects a plurality of portions of the second beam of X-rays that pass through 
15 the volume of tissue from the second beam of X-rays. The X-ray detector generates a second X- 
ray image signal from the plurality of portions of the second X-ray beam detected, which is 
preferably a digital X-ray image signal. 

A resultant image signal is generated based upon a comparison between the first and 
second X-ray image signals generated. In a preferred embodiment, the resultant image signal is a 
20 difference image signal generated by subtracting each pixel of image data of the first digital X- 
ray image signal from each corresponding pixel of the second digital X-ray image signal. When 
the digital X-ray image signals are rendered as two-dimensional, and/or three-dimensional 
images, the pixel-by-pixel subtraction will correspond to subtractions at the same spatial location 
within the corresponding two-dimensional and/or three-dimensional image. Although the 
25 difference image signal generated is based upon a subtraction of image values in the preferred 
embodiment, the general concept of the invention is not limited to only subtracting image values. 



21 



Attorney Docket No. 41905-01419 
5 For example, a ratio between corresponding values could be taken, or innumerable other 
mathematical operations could be performed, as suitable, such as dividing all values by a 
normalization factor, etc. Furthermore, groups of pixels may first be combined in various ways, 
such as pixel averaging prior to the generating a difference image signal. The image may be 
preprocessed by two-dimensional and/or three-dimensional Fourier transformation and 
1 0 comparing the images in the spatial-frequency domain either with or without prior spatial 
filtering. 

The difference image signal is rendered as an image, preferably rendered as a real-time 
image on a display screen to provide information during thermal therapy. The displayed image of 1 
the difference signal provides information on the size and/or volume of localized regions of 
15 tissue, such as denatured tissue, as well as the rate of temperature changes or temperature- 
induced changes and the corresponding relative temperatures. The user may employ the 
displayed information to determine whether desired regions of tissue are receiving enough, or too 
much, heat, as well as determining the extent of tissue necrosis due to the thermal therapy. The 
displayed image may be utilized by the user to determine whether to continue the thermal 
> 20 therapy and to determine whether the heat application parameters need to be altered. For 

example, the amount of heat applied . may be increased, or decreased, or the location where the 
heat is being applied may be altered. 

If the user determines that the thermal therapy should continue, heat is again applied to 
the tissue. The heat application may be a continuous process in which the successive applications 
25 of heat may be considered to be contiguous time intervals. Alternatively, the heat may be applied 
in separate intervals, separated by periods in which no heat is applied to the tissue by the heat 
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5 application assembly. 

The region of tissue is then illuminated with another beam of X-rays in which portions 
are detected by the X-ray detector after they pass through the tissue. The beams of X-rays 
illuminating the tissue may be discrete illumination periods separated by non-illumination 
periods, or may be contiguous time periods forming a continuous illumination by X-rays 

10 classified in time intervals. The X-rays detected after the third illumination beam of X-rays lead 
to a third digital X-ray image signal generated according to the preferred embodiment of the 
invention. The same reference X-ray signal is retrieved from the data storage unit as was 
previously retrieved from the data storage unit. In this embodiment of the invention, the same 
reference X-ray image signal is used to generate difference image signals a plurality of times. 

15 This reference X-ray signal is thus called a static reference X-ray signal to indicate that it will 
remain unchanged for at least a plurality of difference image signals generated. 

After retrieving the static reference X-ray image signal from the data 20 storage unit and 
generating a second resultant image signal, preferably a difference image signal,. it is then 
rendered on a display screen. In the preferred embodiment, the user has the displayed image 

20 available to update determinations on size and/or volume of the cancerous and/or necrotic tissue 
as well as a temperature map, or map of temperature-induced changes, and temperature change 
map within the tissue area of observation. This process may then be repeated numerous times 
until the surgeon determines that the thermal therapy should be terminated. 

The general concepts of the invention include generating both two-dimensional, and 

25 three-dimensional X-ray image signals (e.g., stereotactic X ray imaging), and generating the 
corresponding resultant X-ray image signals. In a preferred embodiment, three-dimensional X- 
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5 ray image signals are produced. In the case in which three-dimensional image signals are 

produced by geometric triangulation methods, at least two X-ray beams through a given point 
within the illumination tissue must be detected. In order to obtain such data, the X-ray 
illumination source and/or detector may be moved relative to a tissue being illuminated, or there 
may be a plurality of illumination sources and/or X-ray detection elements. 

10 FIGURE 6 illustrates another preferred embodiment of a method of 10 thermally 

inducing and monitoring changes to localized regions of tissue in accordance with the invention. 
The description of the method of the preferred embodiment in FIGURE 5 carries over to many of 
the steps illustrated in FIGURE 6. In the embodiment of FIGURE 6, the reference signal 
> generated and stored in the data storage unit is a dynamic reference signal which is updated 

15 subsequent to each resultant image signal that is generated. After generating the resultant image 
signal, rendering and displaying the resultant image, determining the physical properties of the 
localized regions of tissue, and deciding to continue with thermal therapy, the second X-ray 
image signal replaces the first X-ray image signal in the data storage unit to become the up-dated 
dynamic reference signal. Heat is then applied to the tissue, followed by illumination of the 

20 - tissue with a third beam of X-rays. The third beam of X-rays is detected after passing through the 
tissue illuminated with the X-rays, and a third X-ray signal is generated. The dynamic reference 
. X-ray image signal is retrieved from the data storage unit and subtracted, pixel by pixel, from the 
third X-ray image signal to generate a second difference image signal. The second difference 
image signal is rendered and displayed on a display unit for a user to determine the various 

25 physical properties of the localized regions of tissue. If the surgeon determines that thermal 

therapy should continue, the third X-ray image signal replaces the previously X-ray signal in the 
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data storage unit to become the up-dated dynamic reference image signal. This process is 
continued repeatedly until the surgeon determines that thermal therapy should be terminated. 

The preferred embodiments illustrated in FIGURES 5 and 6 describe using either a static 
reference X-ray image signal, or a dynamic reference X-ray image signal. The general concept of 
the invention includes a combination of both static and reference image signals in which the 
reference signal may be static for a period of time, dynamically updated and then could become 
static again to provide an unlimited range of possible combinations between the two limits. 

The following subsections describe several experiments conducted with apparatuses and 
methods according to the concepts and embodiments of this invention. 

MAMMO VISION EXPOSURE EXPERIMENTS ON LASER 
TREATMENT OF PORCINE TISSUE 

This subsection describes preliminary imaging experiments using the MAMMOVISION 
System (85200G-2, Fischer Imaging Co. Denver, CO) during laser-irradiation of porcine fatty 
tissue with a 980 nm diode laser (AOC 25, AOC Medical Systems, South Plainfield, NJ). 

The Experimental Method: 

A 20 x 8 x 3 cm piece of fatty porcine tissue with the skin was compressed at room 
■temperature of 25°C with the MAMMOVISION compression paddle. In a first experiment, an 
optical fiber with a diffusing tip with a diameter of 1 .7 mm (REM series B) was inserted in the 
tissue after making a path with a stainless steel trocar. After the initial scout and stereo images 
were taken to ensure that the fiber was within view of the camera, we began the auto-firing 
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sequence. The sequence was set to take exposures every 30 seconds, although the system 
actually took exposures every 32 seconds. The generator settings were 25 kV, 100 mA, and 90 
mAs. The laser was fired continuously at 10 W starting after the first auto-firing exposure for 
approximately 9 minutes. After approximately 2 minutes, we moved the fiber back so that the 
tip would appear in the images. After the 18th exposure we turned off the laser and continued to 
take exposures for 2 minutes. 

In a second experiment a bare fiber with a 0.6 mm diameter was used (3M series B). The 
experimental procedure was the same as with the diffusing fiber. The generator settings were 25 
kV, 100 mA, and 50 mAs. The laser was fired from between the third and 15th exposures and 
again took images for two minutes of the cooling down period. It should also be noted that we 
varied the laser power from 5 to 10 W after one minute and from 10 to 15 W after another two 
minutes. 

In order to measure quantitatively the difference between the exposures, the average 
intensity (in arbitrary units corresponding to X-ray density) in a region of interest surrounding 
the optical fibers (-30 x 20 mm for the diffusing tip and -26 x 1 8 mm for the bare fiber) was 
recorded for each image. 

Results and Conclusion: 

The average intensity of the region of interest increases approximately linearly during 
laser treatment of the tissue (FIGURES 7 and 8). The rates at which the average intensity 
increases are 34 units/mm for the diffusing tip experiment and 54 units/mm for the bare fiber 
experiment. After stopping the laser, the average intensities of the regions of interest decrease, 
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although not to their initial levels. Both the average intensity increase during laser treatment and 
the decrease after laser treatment are higher in the bare fiber experiment. 

Note: 

The increases in the average intensity may either be a result of changes in the tissue 
density or to the camera's sensitivity to heat. 

CONTINUATION OF MAMMO VISION EXPOSURE EXPERIMENTS 
ON LASER TREATMENT OF PORCINE TISSUE 

This section describes a continuation of the preliminary imaging experiments with the 
MAMIMOVISION System (85200G-2, Fischer Imaging Co., Denver, CO) during laser 
irradiation of fatty porcine tissue with a 980 nm diode laser (AOC 25, AOC Medical Systems, 
South Plainfield, NJ). 

The Experimental Method: 

A 20 x 8 x 3 cm piece of fatty porcine tissue with the skin was compressed at room 
temperature of 25 °C with the MAMMOVISION compression paddle. An optical fiber (REM 
Series B) with a diffusing tip of diameter 1.7 mm was inserted in the tissue after making a path 
with a stainless 15 steel trocar (2.1 mm diameter). The laser was fired for 6 minutes at 10 W 
after, the first image was taken and turned off for the remainder of the experiment. 

After the initial stereo images were taken to ensure that the fiber was within view of the 
camera, we began the auto-firing sequence. The sequence was set to take images ever 30 seconds 
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5 for the first 30 exposures, although the 20 system actually took exposures every 32 seconds/At 
this point (approximately minutes into the experiment), the auto-firing sequence was set to take 
exposures every 6 seconds. 

The difference between the exposures was measured quantitatively by obtaining the 
average intensity (in arbitrary units corresponding to the optical density) in different regions of 
10 interests above and below the optical fiber 148 (see, FIGURE 9). Regions directly surrounding 
the fiber were not studied because the fiber appears to shift through the images and its 
appearance in a region where it was not initially placed would skew the average intensity. 

Results: 

15 The following graphs in FIGURE 10 illustrate the average intensity increases during the 

firing of the laser (the first six minutes) and decrease and level off during the remainder of the 
experiment. In general, the regions closer to the optical fiber show a more significant increase in 
average intensity during the laser irradiation than the regions further away. 

20 Conclusion: 

Although the fluctuations in average intensities of the regions of interest appear small, the 
increases correspond with the time that the laser is on. This suggests that the increase in average 
intensity is correlated to the laser irradiation (possibly by its thermal effects). Since the regions 
further away from the optical fiber show smaller increases than the closer regions, it appears that 
25 the effect of the laser irradiation on the tissue decreases with distance. The average intensities 
decrease after the laser is stopped and eventually level off This may be caused by the 
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5 temperature in the tissue dropping after irradiation and reaching an equilibrium. 

MAMMOVISION EXPOSURE STABILITY EXPERIMENTS OF THE 
PMMA CALIBRATION PHANTOM 

1 0 This subsection describes stability experiments conducted with the MAMMOVISION 

System (85200G-2, Fischer Imaging Co., Denver, CO) on the PMMA calibration phantom (10 x 
10 x4 cm) studied to determine if there are changes in average intensity of the 50 x 50 mm 
images produced by the CCD camera with time during the auto-firing sequence. The effect of an 
additional PMMA plate behind the calibration phantom is also presented. 1 

15 

The Experimental Method: 

The 10 x 10 x 4 cm PMMA calibration block at room temperature (25°C) was held in 
place before the MAMMOTEST camera with the compression paddle. In a first experiment, the 
auto-firing sequence was set to 5 take exposures every 30 seconds for approximately 15 minutes 

20 (after the initial scout and stereo images were taken). The initial scout was taken in the 

Autoexposure mode and its parameters were used in the sequential firing. The generator was set 
manually at 25 kV, 280 mAs, and 100 mA for the sequential firing. The first image was 
accidentally taken at automatic generator settings and was not included in the data discussed 
below. In the second experiment, the auto-firing sequence was set to take exposures every 6 

25 seconds for 3 minutes. The generator settings were set manually at 26 kV, 280 mAs, 100 
mA. 

The minimum, maximum, and average pixel value of these images were recorded. The 
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5 maximum pixel intensities of the center 40 x 40 mm area of the 50 x 50 mm whole image was 
then compared to that of the entire image to avoid potential edge effects (see, FIGURE 1 1). A 
Turbo C program was written to determine the value and location of the pixel with the maximum 
intensity. 

The final experiment compared the effect of placing an additional 0.5 cm PMMA plate, 
10 which will be used in the tissue phantom holder, between the camera and the phantom on the 
average pixel intensity. Three images (an initial scout and two stereo images) were taken at 
generator settings of 26 kV, 380 mAs, 100 mA. 

Results: 

1 5 The average intensity of the images in the first and second experiments are within 0.5% 

and 0.2%, respectively (FIGURES 12 and 13). The standard deviations of the average intensities 
are too small to be visible in FIGURES 12 AND 13. The maximum pixel intensities in both 
experiments show large variations. Much of these variations disappear when analyzing the center 
40 x 40 mm region, suggesting that the artifact causing these variations occurs more frequently 

20 in the edges of the images (see, FIGURES 14 and 15). In the five images where the variations 

were still apparent in the 40 x 40 mm region, the variations do not appear in a 30 x 30 mm center , 
region. The minimum intensity values remained virtually unchanged in the center region (Tables : 
1-2). According to Tables 3-4, the maximum intensity occurs on the same few pixels to 
the experiment. In both experiments, the maximum intensity value occurs mostly on pixel : 

25 coordinates (628, 1 020, maximum intensity 3 1 70) and ( 1 0 1 8, 1 79, maximum intensity 2400). 

The average of the average intensities of the images taken in the last experiment with the 
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additional PMMA plate is 1393, or 66% of the value without the plate, which indicates that the 
plate introduces a significant attenuation. However, this effect may be compensated for by using 
other X-ray generator parameters. 

Conclusion: 

These experiments demonstrate that the pixel intensities of the sequential images are 
stable, with the exception of a few pixels. 

Note: 

The minimum pixel intensity is zero and corresponds to total attenuation (i.e., no X-ray 
photons reaching the camera) and a white image. The maximum pixel intensity is 4095 and 
corresponds to no attenuation (i.e., all X-ray photons reaching the camera) and a black image. 
Each pixel corresponds to an area of approximately 50 x 50 pm. 
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Table 1 : Pixel intensity data of Experiment 1 (30 s firing rate, 25 kV, 280 mAs, 100 mA) 

(RMS [root mean square] represents the standard deviation of the intensity 
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Table 2: Pixel intensity data of Experiment 2 (6 s firing rate, 26 kV, 280 mAs, 100 mA) 

(RMS [root mean square] represents the standard deviation of the intensity 

over the whole image) 
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Table 3: Values and pixel coordinates of maximum intensities 
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ANALYSIS OF THERMALLY-INDUCED CHANGES 
IN MAMMOVISION IMAGES BY PIXEL A VERAGING 

This subsection describes test results that display changes in sequential images taken with the 
MAMMOVISION System on the PMMA calibration phantom and on laser-irradiated fatty 
porcine tissue. 

The images taken with the MAMMOVISION system are 50 x 50 mm images containing 
1024 x 1024 pixels with each pixel corresponding to a 50 x 50 p.m area. Since the desired target 
volume is 30 mm in diameter, the main region of interest (ROl) during treatment is a 30 x 30 
mm area in the image. These data in ASCII (text) files containing the intensity for each pixel are 
5 MB each. These files are too large to handle when dealing with multiple images. In addition, 
the resolution of the images is much greater for each image than is necessary for thermal 
imaging; where a resolution on the order of 1 x 1 mm is sufficient. It is therefore beneficial to 
convert the images from 1024 x 1024 pixels to a smaller resolution by averaging blocks of the 
1024 x 1024 into one pixel (see, Table 4). Another benefit of reducing the resolution is 
removing the noise effects generated by a few unstable pixels. 
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5 Table 4:Number of pixels and resolution of averaged images 



Number of Pixels 


Size Averaged 
Block 


Pixel Size of 
Average Images 


Number of Pixels in 
30 x 30 mm ROI 


1024 x 1024 
(no averaging) 


1 x 1 


50 x 50 /jm 


600 x 600 


512x512 


2x2 


0.1 x 0.1 mm 


300x300 


256 x 256 


4x4 


0.2 x 0.2 mm 


150 x 150 


128 x 128 


8x8 


0.4 x 0.4 mm 


75x75 


64x64 


16x 16 


0.8 x 0.8 mm 


37x37 



Stability of 1024 x 1024 Images: 

In a previous experiment (see, above), the average pixel intensity of the entire 1024 x 
1024 image of the PMMA calibration phantom is shown to be stable within 0.5% in the 
sequential-firing mode. The purpose of this study was to quantify the stability of each pixel. 

The Experimental Method: 

The PMMA calibration block (10x10x4 cm) was held in place before the 
MAMMOVISION (85200G-2, Fischer Imaging Co., Denver, CO) camera with the compression 
paddle at room temperature (25 °C). The initial scout was taken at 00 in the Autoexposure mode 
and its parameters were used in the subsequent firing. The generator was set manually at 25 kV, 
280 mAs, and 100 mA for the sequential firing. The auto-firing sequence was set to take 
exposures every 30 seconds for approximately 14 minutes for a total of 27 images. A Turbo C 
program was written to find pixels that changed in intensity by more than 20% between two 
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successive images. 

5 

Results: 

The analysis of the 27 1024 x 1024 images revealed that only 90 pixels 15 showed 
significant variation (>20%) at least once during the experiment (Table 5). The maximum 
number of pixels that changed by more than 20% from one image to the next was 19. This 
10 analysis corroborates the earlier findings that the pixel intensities of the sequential images are 
stable, with the exception of a small percentage of pixels. 



Table 5: Coordinates of pixels with intensity change greater than 20% 



X- 

coordinate 


Y- 
coordinate 


Number Of 
Occurrences 


Maximum 
Percent 
Change 




X- 

coordinate 


Y- 

coordinate 


Number Of 
Occurrences 


Maximum 
Percent 
Change 




538 


3 


23.2 




711 


1011 


4 


131.5 


157 


110 


2 


22.6 




716 


12 


1 


23.4 


162 


534 


10 


38.8 




731 


805 


2 


34.4 


181 


473 


4 


44.5 




743 


241 


2 


22.2 


208 


521 


10 


30.1 




745 


646 


8 


46.5 


212 


584 


14 


32.6 




749 


358' 


5 


40.0 


227 


23 


3 


25.1 




774 


970 


2 


23.1 
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X- 

coordinate 


Y- 

coordinate 


Number Of 
Occurrences 


Maximum 
Percent 
Change 




X- 

coordinate 


Y- 
coordinate 


Number Of 
Occurrences 


Maximum 
Percent 
Change 


234 


992 


1 


-30.4 




778 


703 


2 


22.2 


253 


516 


1 


20.4 




789 


938 


6 


27.3 


302 


275 


6 


-39.6 




808 


527 


1 


44.2 


305 


652 


4 


23.1 




809 


415 


3 


22.8 


312 


601 


7 


52.3 




825 


586 


2 


24.9 


327 


474 


6 


40.9 




830 


238 


1 


20.2 


329 


1004 


6 


34.6 




847 


127 


1 


21.3 


340 


768 


3 


50.0 




854 


630 


8 


44.8 



X- 

coordinate 


Y- 

coordinate 


Number Of 
Occurrences 


Maximum 
Percent 
Change 




X- 

coordinate 


Y- 

coordinate 


Number Of 
Occurrences 


Maximum 
Percent 
Change 


344 


408 


1 


24.6 




892 


604 


4 


32.2 


389 


885 


3 


75.5 




903 


670 


1 


20.7 


395 


4 


7 


36.2 




903 


832 


4 


26.7 


396 


195 


2 


115.6 




909 


714 


1 


21.5 


421 


882 


1 


. 22.8 




912 


208 


1 


20.5 


457 


899 


5 


28.9 




924 


663 


8 


92.5 


459 


906 


3 


38.3 




926 


164 


5 


48.4 


472 


488 


4 


33.6 




928 


468 


3 


22.7 
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X- 

coordinate 


Y- 

coordinate 


Number Of 
Occurrences 


Maximum 
Percent 
Change 




X- 

coordinate 


Y- 

coordinate 


Number Of 
Occurrences 


Maximum 
Percent 
Change 


476 


324 


1 


-20.8 




929 


983 


7 


75.4 


484 


304 


2 


24.2 




932 


474 


6 


38.5 


484 


584 


1 


20.4 




935 


986 


3 


28.3 


493 


178 


1 


25.8 




942 


262 


3 • 


28.8 


494 


770 


4 


66.6 




947 


400 


2 


20.8 


528 


113 


1 


-21.7 




964 


65 


6 


59.7 


536 


853 


4 


61.4 




966 


681 


1 


20.1 


571 


217 


2 


56.2 




970 


63 


5 


89.5 



X- 

coordinate 


Y- 

coordinate 


Number Of 
Occurrences 


Maximum 
Percent 
Change 




X- 

coordinate 


Y- 
coordinate 


Number Of 
Occurrences 


Maximum 
Percent 
Change 


572 • 


1001 


2 


22.5 




976 


958 


3 • 


26.9 


603 


162 


2 


31.0 




978 


25 


1 


20.7 


614 


906 


6 


32.4 




979 


161 


1 


27.9 


624 


153 


1 


20.2 




981 


234 


1 


27.2 


628 


1020 


10 


99.6 




986 


281 


1 


24.7 


635 


724 


1 


23.8 




987 


. 878 


2 


21.3 


653 


66 


4 


25.5 




994 


221 


8 


42.1 


660 


978 


4 


35.3 




995 


53 


2 


-28.6 
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X- 
coordinate 


Y- 

coordinate 


Number Of 
Occurrences 


Maximum 
Percent 
Change 




X- 

coordinate 


Y- 

coordinate 


Number Of 
Occurrences 


Maximum 
Percent 
Change 


669 


303 


2 


34.4 




997 


91 


3 


23.2 


669 


304 


2 


123.3 


• 


100 


791 


1 


24.5 


669 


388 


2 


22.3 




1003 


809 


7 


36.5 


670 


304 


2 


62.2 




1014 


669 


4 


25.6 


674 


388 


1 


23.7 




1017 


280 


1 


27.7 


682 


0 


4 


83.9 




1018 


179 


2 


56.1 


686 


947 


2 


22.9 




1022 


898 


13 


31.0 



Table 5: Coordinates of pixels with intensity change greater than 20% 

Stability of 128 x 128 Images 

This subsection describes results of quantifying the stability of the 5 pixel intensities of 
128 x 128 averaged images. 

The Experimental Method: 

The images used in the previous analysis were converted from 1024 x 1024 pixel images 
into 128 x 128 pixel images by averaging 8x8 pixel blocks into one pixel, by using a Turbo C 
program. These images were then analyzed 10 by another Turbo C program to find pixels that 
changed in intensity by more than 2%. 
5 Results: 
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5 In the 27 128 x 128 pixel images, only 8 pixels showed variation of greater than 2% (see, 

Table 6). The maximum variation was 3.3%, which was 15 in a pixel that increased by that 

percentage in one image and then decreased to approximately its original value in the next 

image. No other pixel showed more than a 2% change more than once. The maximum number 

of pixels changing by more than 2% between two consecutive images was 2. These results 

10 indicate that any image changes of more than 5% in 128 x 128 images 20 will be significant 

(above noise level) during laser-irradiation. 

Table 6: Coordinates and image number of pixels with over 2% 
intensity change in the 128 x 128 images 



X-coordinate 


Y-coordinate 


Image Number 


Percent change 


6 


73 


4 


-2.04 


10 


59 


15 


-2.07 


35 


96 


20 


2.13 


49 


24 


15 


2.16 


78 


127 


24 


2.15 


83 


38 


11 


-3.12 


83 


38 


12 


3.44 


92 


5 


3 


2.17 


120 


73 


17 


2.03 



15 

Colorized Display of Changes in Laser-Irradiated Tissue: 

This subsection describes displaying the changes in X-ray intensity of laser-irradiated 



44 



4 

Attorney Docket No. 41905-01419 

5 tissue in color according to this invention. 

The Analytical Method: 

The images of the laser-irradiated fatty porcine tissue were obtained from a previous 
experiment. A 20 x 8 x 3 cm piece of fatty porcine tissue with the skin was compressed at room 
10 temperature of 25 °C with the MAMMOVISION compression paddle. An optical fiber (REM 
Series B) with a diffusing tip of diameter 1.7 mm was inserted in the tissue after making a path 
with a stainless steel trocar (2.1 mm diameter). The laser was fired for 6 minutes at 10 W after 
the first image was taken and turned off for the remainder of the experiment. 

After the initial stereo images were taken to ensure that the fiber was 5 within view of the 
15 camera, we began the auto-firing sequence. The sequence was set to take images ever 30 seconds 
for the first 30 exposures, although the system actually took exposures every 32 seconds. 

These images were then converted to 128 x 128 pixel images and each image was 
subtracted from the first image. These images were then colorized using MICROCAL ORIGIN 
to show changes around the fiber with time. Two color scales are compared: linear and 
20 logarithmic. 

Results: 

In FIGURES 17A-21B, the changes in the tissue surrounding the fiber are evident. The 
changes are more apparent in the logarithmic scale than in the linear scale. Movement of the 
25 fiber during the experiment cause the red and white regions at the original location of the fiber 
and the dark blue region at the actual location of the fiber. 
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5 

Conclusion 

From the stability analysis of the 128 x 128 pixel images of the calibration phantom, it 
can be assumed that intensity changes above 5% are assumed to be due either to thermal 
expansion or to tissue movement. These changes show that X-ray images can be utilized to 

10 monitor tissue changes during laser irradiation in accordance with this invention. 

Although the images presented here are in a logarithmic scale to 25 highlight the changes 
around the fiber, the scale will not amplify the temperature range of most interest to 
hyperthermia when the intensity changes are correlated with temperature (see, FIGURES 16A, 
16B and 16C). The temperature range that needs to be monitored most closely in this 

15 embodiment of the invention is between 55 and 65 °C, which is expected to be in the middle or 
higher (red) end of the range of intensity changes. Temperatures above 65°C do not need to be 
monitored by temperature. In order to amplify the correct temperature range, a more suitable 
transfer function may be developed to be converted to a color scale (see, FIGURES 16A-21B). 

20 Note: 

The minimum pixel intensity is zero and corresponds to total attenuation (i.e., no X-ray 
photons reaching the camera) and a white image. The maximum pixel intensity is 4095 and 
corresponds to no attenuation (i.e., all X-ray photons reaching the camera) and a black image. 
After the images are subtracted, the maximum possible intensity range is from -4095 to 4095. 
25 Each pixel corresponds to an area of approximately 50 x 50 /u m in the 1024 x 1024 pixel images 
and to approximately 0.4 x 0.4 mm in the 128 x 128 pixel images. 
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STABILITY TEST OF M AMMO VISION EXPOSURES OF PORCINE TISSUE 
This subsection presents an analysis of a stability test of the MAMMOVISION System 
(85200G-2, Fischer Imaging Co., Denver, CO) exposures on fatty porcine tissue during the auto- 
firing sequence in accordance with this invention. 

The Experimental Method: 

An approximately 20 x 8 x 3 cm piece of fatty porcine tissue with the 25 skin was 
compressed at room temperature 20°C with the MAMMOVISION compression paddle. The 
initial scout was taken in the Autoexposure mode and its parameters were used in the sequential 
firing. The generator was set manually at 25 kV, 50 mAs, and 100 mA for the following images 
in the sequential firing. The auto-firing sequence was set to take exposures every 30 seconds for 
approximately 15 minutes (after the initial scout and stereo images were taken). 

The images obtained from this sequence were converted from 1024 x 1024 pixel images 
into 128 x 128 pixel images by averaging 8x8 pixel blocks into one pixel, by using a Turbo C 
program. These images were then analyzed by another Turbo C program to find pixels that 
changed in intensity by more than 2%. 
Results: 

In the 32 128 x 128 pixel images, only 3 pixels showed variation of greater than 2% (see, 
Table 7). The maximum variation was 8.8%. All three pixels increased by more than 2% in one 
image and then returned to approximately the original value in the following image. These 
results indicate that any image changes of more than 5% in the 128 x 128 images will be 
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significant (above noise level). 

Table 7: Coordinates and image number of pixels with over 2% intensity 
change in the 128 x 128 images 



X-coordinate 


Y-coordinate 


Image number 


Percent change 


111 


90 


6 


8.79 


111 


90 


7 


-8.16 


15 


70 


8 


5.00 


15 


70 


9 


-4.95 


107 


18 


24 


4.11 




X-coordinate 


Y-coordinate 


Image number 


Percent change 


107 


18 


25 


-3.98 



Conclusion: 

This analysis corroborates the earlier findings that the 128 x 128 pixel images are stable, 
although those findings indicated that there were no changes above 5% while the results from 
this analysis indicate that one pixel changed above 5%. This discrepancy may be explained by 
the relatively low generator mAs value (50 mAs) in this experiment compared to that of the 
previous experiment (280 mAs), since increasing the mAs value increases contrast. Changes 
above 5% can be assumed to be due either to thermal expansion or to tissue movement. 

Note: 

The minimum pixel intensity is zero and corresponds to total attenuation (i.e., no X-ray 
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photons reaching the camera) and a white image. The maximum pixel intensity is 4095 and 
corresponds to no attenuation (i.e., all X-ray photons reaching the camera) and a black image. 
Each pixel corresponds to an area of approximately 50 x 50 ju m in the 1024 x 1024 pixel images 
and to approximately 0.4 x 0.4 mm in the 128 x 128 pixel images. 

Although only the presently preferred embodiments have been described in detail above, 
people with ordinary skill in the art will readily appreciate from the teachings herein that many 
modifications are possible in the exemplary embodiments without materially departing from the 
novel teachings and advantages of this invention. 
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